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Abstract 



The interfacial profile between coexisting phases of a binary mixture (A,B) in 
a thin film of thickness D and lateral linear dimensions L depends sensitively on 
both linear dimensions and on the nature of boundary conditions and statistical 
ensembles applied. These phenomena generic for systems in confined geometry 
are demonstrated by Monte-Carlo simulations of the bond fluctuation model of 
symmetric polymer mixtures, using chains containing Na = Nb = N = 32 effec- 
tive monomers connected by effective bonds with an attractive interaction between 
monomers of the same type and a repulsive interaction between different types. We 
use short range potentials at the walls, the right wall favoring A monomers and the 
left wall B monomers. Periodic boundary conditions are applied in the directions 
parallel to the walls. Both the canonical and semi-grand-canonical ensemble are 
studied. We argue that the latter case is appropriate for experiments with a lateral 
resolution L much less than the actual lateral sample size, in thermal equilibrium. 
In the canonical ensemble, the interfacial width w increases (from small values which 
are of the same order as the "intrinsic profile" ) like w oc \^D, before a crossover to a 
saturation value w max (w^ ai oc In L) sets in. In the semi-grand-canonical ensemble, 
however, one finds the same widths w oc \[~D as in the canonical ensemble for not 
too large L, while for large L the interfacial profile is smeared out over a finite 
fraction of the film thickness (w oc D for D — > oo). We discuss the implications of 
these findings for the interpretation of both simulations and experiments. 



1 Introduction 



Since van der Waals' famous treatment of the interface between coexisting fluid and 
gas, calculations and measurements of interfacial profiles between coexisting phases have 
remained a particular challenge 0, ||, |3|, ||, ||. Van der Waals type mean field theories 
were extended to mixtures by Cahn and Hilliard ||, as well as to more complex systems 
such as polymer solutions 0, polymer blends || |9|, |10 



microemulsions [O, tL2L etc 



These theories, but also more complex ones such as the self-consistent field approach 
to strongly segregated polymer mixtures |13|, [L4], [15], [L6|, [T7|, [Tj| [19], |20|, ^Tj, |22| consider 
a hypothetical "intrinsic" profile that is neither directly accessible by experiment nor by 
computer simulation. I.e., both in experiment and in simulations these "intrinsic" profiles 
- if they have a well-defined meaning at all — are broadened by fluctuations (such as 
capillary waves \W , |20| , p3| which can be modelled by harmonic displacements of the local 
interface position II, ||, HE), ^ s broadening is also limited by the geometry of 
the (finite !) system that is considered and affected by forces due to the boundaries of 
the system. In fact, the latter may give rise to interface binding/unbinding phenomena 
("wetting" ^ ^B], |27], [HJ), and the associated fluctuations do have a strong 

effect on the observable average interfacial profiles [31 , [H| [33], |34], |35| . (Note that external 



potentials such as gravity, which clearly limit interfacial fluctuations in liquid-gas systems, 
are disregarded here throughout.) 

There are numerous examples in the literature, both in experiment (e.g. on polymer 
8|) and simulation (ranging from lattice gas models to the water-vapor 
where this fluctuation broadening of interfacial profiles is disregarded and 
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the resulting interfacial width is interpreted as if the intrinsic width were obtained. Recent 
however, indicates that such an approach may yield rather misleading results ! 
Fig.l sketches the L x L x D geometries that are typically used in simulations |23|, 
[33| , [34] , [35| , |39| , fdOfl and experiments @, [37], |38|. L is the linear dimension parallel 
to the interface(s), D perpendicular to the interface(s). Fig. la shows a system with a 
single interface between two demixed phases using either hard walls as boundary condi- 
tion normal to the interface (left part) or "antiperiodic" boundary conditions (i.e., an A 
monomer leaving the right boundary of the simulation box reenters at the left boundary 
as a B monomer, and vice versa; right part). This antiperiodic boundary condition is 
particular useful for strictly symmetrical mixtures where two species A,B intrinsically are 
not different apart from their "label" , and thus one can achieve a geometry with a single 
interface and no disturbing effects due to external walls. If periodic boundary conditions 
are used in all three directions, the system must contain two interfaces (Fig. lb). We 
emphasize that the concentration profile (and hence also its width to) obtained in these 
situations depends both on D and on L, and in general may differ appreciably from any 
sensibly defined intrinsic profile ! The present paper intends to study these dependences 
for a simple model system that is nevertheless close to physical systems of practical inter- 
est, namely polymer mixtures |l 



14], [37], pq| . Experiments are done in the geometry of 
Fig.lc and encounter closely related problems of interpretation as the simulations as will 
be explained below. 

While in the thermodynamic limit the different ensembles of statistical mechanics 
yield equivalent results, this is not true in systems which have all their linear dimensions 



finite |41 |. Considering binary (A,B) mixtures at constant density, we may deal either 
with the canonical ensemble where the average volume fraction $ of species A is held 
fixed (if PaiPb are the average densities of A particles/B particles in the system, then 
$a = $ = Pa/(Pa + Ab); an d $b = 1 — $ = Pb/(Pa + Pb) being the volume fraction 
of species B) or with the semi-grand-canonical ensemble where the chemical potential 
difference Ap between A and B is held fixed. Note that in the left part of Fig. la the use 
of the canonical ensemble suppresses the fluctuation of the average position of the interface 
- e.g. for $ = 1/2 and "antisymmetric" walls acting on a "symmetric" mixture |42| (as 
considered here and in Refs. |32|, the average position of the interface is fixed in the 



middle of the film, z = D/2. (We disregard here the case [[31], [32], |33| of temperatures below 
the transition temperature T C (D) of the interface localization-delocalization transition 
where the interface is bound either to the left or the right wall and $ = 1/2 is maintained 
only via a lateral inhomogeneity, i.e., formation of domains separated by interfaces running 



perpendicular to the walls ||43|| .) In contrast, using the semi-grand-canonical ensemble, 
$ and hence the average interface positions are fluctuating and, in particular for small 
L, this fluctuation of the average interface position may even be more important than 
local distortions of the interface. Naturally, one has to be careful in the interpretation 
of experiments, too: While the canonical ensemble is appropriate for a discussion on the 
length scale L samp i e comprising the total sample, this is not true if the measurement is done 
with a technique characterized by a length scale L of lateral resolution with L <C L samp i e 
(Fig.lc). The condition $ = const fixes the average interface position on the length scale 
Lsampie, but it clearly does not fix it on the length scale L that is probed. This is the case 



for experiments where interfacial properties are measured "in situ" in the melt phase |44 
We may consider the slab of width L that is measured as a "subsystem" in a statistical 
sense that can exchange A,B particles freely by diffusion with the rest of the system. This 
acts as a reservoir, in thermal equilibrium, and hence the semi-grand-canonical ensemble 
applies. Implicit in this interpretation is the assumption that the time constant of the 
measurement is large enough such that one measures not just a kind of "snapshot" of 
the interfacial configuration, but carries out a meaningful statistical average in thermal 
equilibrium. 



For the fully periodic system [35, £|0) (Fig. lb), on the other hand, semi-grand-canonical 
techniques are inappropriate since then the two-phase configuration shown is at best meta- 
stable: If one simulated long enough, the two interfaces would meet and annihilate each 
other, and either the A-rich domain or the B-rich domain would vanish altogether. The 
condition $ = const in the canonical ensemble and fixing the center of mass of the B- 
rich phase then suppresses also fluctuations of the average interface positions. Finally, in 
the antiperiodic case there is a clear translational invariance; no origin at the z-axis is 
distinguished except by the location of the interface itself. Then it is natural to record 
all information in a moving frame: In each configuration that is analyzed the point z = 



is chosen to coincide with the interface position averaged in x,y directions [£3|, |45[j. Even 
in the semi-grand-canonical ensemble, the fluctuation of the average interface position is 
thus eliminated per construction. 

In the present paper we shall restrict our attention to the case of a single interface 
confined between "antisymmetric" walls (i.e., the left part of Fig. la) since this situation 
is closest to the experimental one (Fig.lc), studying the finite size effects both in the 



canonical and the semi-grand-canonical case. We anticipate, however, that many of our 
findings carry over qualitatively to the other cases (antiperiodic boundary conditions in 
Fig. la, periodic boundary conditions in Fig. lb) as well. 

In the next section, we define the model that is used in our study and make a few 
remarks on the simulation technique. Sec.|3] briefly reviews the theoretical background, 
while Sec.fl gives our results on interfacial profiles and interfacial widths and discusses the 
correlation function of interfacial fluctuations in directions parallel to the walls and the 
associated correlation length £||. Sec.|5] describes a test of capillary wave concepts. Finally, 
Sec.^j summarizes our conclusions. 



2 Model and simulation technique 

Since we desire to make as close contact as possible with experiments on polymer mix- 
tures confined in thin film geometry, where size effects as discussed here have already 
been observed |34|, [|7[], we study a coarse-grained model for a symmetric polymer mix- 
ture, extending previous studies of the bond fluctuation model on the simple cubic 
lattice |E| |23], |50|, |51}| . We summarize here only very briefly some key features |52| 
of this well studied model: One integrates n = 3 — 5 chemical bonds along the back- 
bone of a real polymer chain into one effective bond, whose length is allowed to fluctuate 
between the following values: 2, a/5, a/6, 3, vTU, all lengths being measured in units of 
the lattice spacing. This variation in the bond length is not only computationally ad- 
vantageous |52|, |53| , but also physically reasonable as a result of the internal degrees 
of freedom of the combined 3 — 5 original chemical bonds. Furthermore, the 87 differ- 
ent bond angles between subsequent effective bonds approximate properties of continuum 
models already rather closely. These effective bonds are joined at effective monomers, 
each blocking all 8 sites of an elementary cube from further occupation. We work at a 
volume fraction = 1/2 of occupied lattice sites, which corresponds to a dense melt ||53|| , 
i.e., the density of effective monomers is pb = 1/16 in our units. For consistency with 



previous work [EJ where interfaces for the antiperiodic boundary conditions (Fig. la, right 
part) were studied we choose here a chain length N = 32, which hence corresponds to a 
degree of polymerization N p = nN « 96 — 160. For this choice, the radius of gyration 
R g of the chains is about R g ~ 7 lattice spacings. Defining then the effective statistical 

segment length b from the relation valid for Gaussian chains, R g = b^jN/6, we obtain 
b « 3.05. Note that the end-to-end distance is R ~ 17, i.e., within our accuracy the rela- 
tion R = b\^N is also fulfilled. Of course, it would be nice to study much longer chains so 
that the length scales b and R g are more distinct from each other, but this would require 
an unreasonable large effort in computer time and thus has not been attempted in the 
present context. 



The Hamiltonian is written in an Ising model type spin representation H2 



as 
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where Si — +1 if monomer i is of type A, Si - 
parameter e being defined as k B Te = eaa = £bb 



— 1 if it is of type B, the interaction 
= — €ab, if the distance d(i,j) between 



monomers does not exceed a/6, while for more distant neighbors the interaction is zero. 
From previous work |50| it is known that the critical point of phase separations occurs at 
tcrit = 0.0144, and since we wish to avoid critical fluctuations, we work here at e = 0.03, 
i.e., far below the critical point (T/T crit = 0.48), although experiments are done for [34] 
T/T crit pa 0.95. In the latter case, however, the correlation length £ of order parameter 
fluctuations is already about twice as large as the radius of gyration, while in our case it is 
distinctly smaller, £ pa 3.6. Since the correlation length sets the scale for the concentration 
profile — at least in its wings || — one would need to simulate much larger systems for 
T/T crit pa 0.95 than for the present choice, and in fact with the present computer facilities 
it is not yet possible to work so close to the critical point. Our choice corresponds to 
a Flory-Huggins parameter [55, [5I| \ = ^ z eff e ~ 0.159 since the effective coordination 



number is |23| z e ff ~ 2.65, which is still physically reasonable although it corresponds to 



a strongly segregated situation (%iV pa 5.09). From previous work the interfacial tension 
E is known |23[] as well; namely E//c B T = a ~ 0.015. 

In Eq.|l], a wall interaction of square-well type has been added if monomers are not 
more distant from the wall than 2 lattice spacings. We choose e w = 0.1 since from 
studies of wetting phenomena in polymer mixtures ]57 we can safely expect that this 
choice leads to a wetting transition temperature T w below the considered temperature. 
Note that the transition temperature T C (D) of the interface localization-delocalization 



transition occurs 0, |33|] is expected close to T w , and we need T > T C (D) in order to 
have the situation sketched in Fig. la (interface delocalized and on average in the middle 
of the thin film). The absolute magnitude of both surface terms is chosen the same, the 
sign opposite, so we realize an antisymmetric wall situation, as desired. 

The linear dimensions L parallel to the walls are in the range from L = 64 to L = 512, 
while 8 < D < 64. Remembering that R g ~ 7, our film thicknesses thus range from one 
R g to about 9R g . 

Since these system sizes are rather large, at least for the higher values of D and 
L, and interfacial fluctuations are notorious slow p2 |, one must pay much attention to 
an adequate equilibration of the model system. In the canonical ensemble the order 
parameter is conserved and we have ~p A = ~p B . Two different types of runs have been 
made: In the first type of runs, "slithering snake" moves |52| (a chain end of a randomly 
chosen chain is cut at one end and one attempts to join it at the other end in a randomly 
chosen direction) are mixed with the "random hopping" moves ||52| , [5311 (one tries to move 
a randomly chosen monomer of a randomly chosen chain in a randomly chosen lattice 
direction by one unit). It is clear, however, that long range concentration fluctuations 
will relax rather slowly with this algorithm — diffusion of the chains through the whole 
system is required. To ease this problem, a second type of runs was tried where the 
two types of moves mentioned above were combined with chain identity exchanges. To 
conserve the order parameter, one needs to transform a randomly chosen chain (1) from 
A to B and another randomly chosen chain (2) from B to A. After some tests we found 
it useful to mix these moves (slithering snake : random hopping : identity exchange) in 
the ratios 3:1:0.1 (note that the "1" for the random hopping means one attempted move 
per monomer of the chosen chain). The total amount of these moves is then counted 
as "4 Monte-Carlo steps (MCS)" in order to define a "pseudotime" unit for this mixed 
algorithm. For the semi-grand-canonical ensemble, we always use an algorithm of this 



second type, the only distinction being that one attempts an identity change of a single 
chain (A to B or B to A, respectively) instead of an identity exchange between two chains. 
Since the volume fraction is not conserved, p A = p B holds on average only, but in general 
not for individual configurations generated. 

Of course, a price one has to pay for these mixing of algorithms is that one performs 
considerably less moves per second than with algorithms that carry out reptation moves 
only. E.g., on an IBM RISC 6000-250 workstation one needs only 30 sec to perform 10 6 
MCS per chain if only reptation moves are carried out, while the canonical moves including 
chain exchange take t ~ 140 sec, the semi-grand-canonical moves even t pa 150 sec. In 
order to run a system with L = 512 and D = 32 for 10 6 MCS for every chain in the system 
with semi-grand-canonical moves, 700 hours CPU are needed. Since we found it necessary 
to carry out runs that were partly up to 10 times longer, and many combinations of (L, D) 
need to be studied, this project is not feasible at workstations as quoted. Fortunately, 
more powerful machines (including DEC 8400 and SG PowerChallenge) could be accessed. 



3 Theoretical background 
3.1 Definitions 

In this section we define the quantities that will be calculated in the later sections and 
recall some of the basic theoretical predictions. We denote the coordinate of a point inside 
two walls as (r, z) where r = (x, y) is a two-dimensional coordinate in the surface plane 
of the left wall (at z = 0), and z measures the distance from this wall. From the local 
densities pa{t, z), Pb{t, z) of all monomers we define then a local order parameter m(r, z) 
as follows: 

M (v)- riv !' f, ! f 1 , (2) 

p A (r,z)+p B (r,z)' 
while the average order parameter profile is 

m(z) = L~ 2 / dr m(r, z). (3) 



Here / dr stands symbolically for a summation over the L 2 lattice points in plane z, and 
m(z) is defined at the discrete points z = 1, 2, D — 1 only, of course. 

A quantity of key interest is the order parameter correlation function in the direction 
parallel to the walls, 

= (m(r,z) m(0,z)) - (m(0,z)) 2 
9[ ' ' (m(0, zf) - (m(0, z)) 2 ' 1 1 

from which we can extract a correlation length £|| by observing at large distances r an 
exponential decay, 

g(r,z = D/2) oc -j^exp (-^j ■ ( 5 ) 

Here we have anticipated (as is also borne out by our numerical data presented below) 
that the slowest decay of g(r, z) occurs right in the center of the film, due to the interfacial 
fluctuations |1], [3|, 0. 



We define the interfacial width w from a fit to a tanh profile, 
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where the densities of A monomers at the coexistence curve p^Acoexi P^Acoex are known 
from previous work 50[| . Alternatively, we can write Eq.^ in terms of m(z) in a more 
compact form 
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(7) 



+ Pa coed > noting that for our symmetric mixture 

Actually, for our choice of e = 0.03 the mixture in 

the bulk is essentially fully segregated [^3| 
nib ~ 1, and hence Eq.^| simplifies to 



n {1) 
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Ph,coex ~ 



1 + tanh 



Since we have chosen to work with only even values for D, the "Gibbs dividing surface" ]T|, 
U at z = £)/2 always coincides with a position at the lattice. 

The fact that the bulk order parameter m& ~ 1 in our model implies that in the 
present case the bulk order parameter fluctuations are negligible, and any interplay of 
bulk and interfacial fluctuations also can be disregarded here. In this sense, the present 
model is simpler than the Ising model considered in Refs.[32, 33 1. (One cannot go to 
such low temperatures with nib ~ 1 there because the interface roughening transition [ 58 ] 
intervenes.) 



3.2 Capillary waves 



Our choice of the tanh profile in Eqs.|6|-|8J is motivated by the fact that this form is predicted 
by mean field type theories, such as Cahn-Hilliard type long wavelength theories [|7|, |pj 
and the self-consistent field theory in the limit [13, 14[ xN ~^ 00 • All these theories, 
however, neglect fluctuations in the position of the "local" Gibbs dividing surface which 
we may define as [pa,b(x, y) = D" 1 fdz pa,b(%, U, z)\ 
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Pa(x,v) + Pb(x,v) 
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It then is convenient to consider the deviation h(x, y) = z int (x, y) — D/2 of this fluctuating 
interface from its average location. As is well known 0, |3|, |J, the idea of capillary wave 
theory is to put the free energy cost of these fluctuations proportional to the increase in in- 
terfacial area caused by these fluctuations. Assuming that the gradients \dh/dx\, \dh/dy\ 
are small, this yields (normalizing again by a factor k B T) 
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(10) 



where it is assumed that a is simply the bulk interfacial tension (normalized per tempera- 
ture). By Fourier transformation one obtains from a simple Gaussian Hamiltonian, 
h(q) denoting the Fourier component of h(x, y) for wavevector q, 

^F cw = -Y j q 2 \Hq)\\ (11) 

£ q 

and from the equipartition theorem one can immediately conclude that the mean square 
value of h(q) is @ 

(\Hq)\ 2 } = — 2 . (12) 
aq z 

This power law spectrum clearly leads to divergences when one computes the local mean 
square displacement of the interface jl], 0, 0, ffl, 

s 2 = (h 2 (x,y)) 

= E(l%)l 2 > 



dq(\h(q)\ 2 ) 




(13) 

Since the integral fdq/q diverges logarithmically both for q — > and q — > oo, we have 
heuristically used both a lower cut-off q min and an upper cut-off q max . It is clear that 
the naive approximation leaving the 1/q 2 spectrum in Eq.[12] completely unmodified in 
between q min and q max and using these sharp cut-offs needs a closer investigation. It could 
be necessary to cut off these divergences in a smooth way, using suitable correction terms 
in the capillary wave Hamiltonian, EqsjH] and 11, in order to obtain a more accurate 



description |p9|| . Even if one accepts Eq.[L3], the correct choice of q min and q max is a 
problem, particularly for a polymer mixture. For large enough D such that no confinement 
effects are felt by these interfacial fluctuations, the minimal value of q possible in the 
geometries of Fig.la,b clearly is q m i n = 2n/L (remember that we work with periodic 
boundary conditions), but the choice of q max is much more questionable: For a small 
molecule system away from the critical region, one just takes q max = 2n/a, a being a 
molecular diameter, while near the critical point q m ax = 27r/£ as the correlation length 
£ then is the only important length scale in the system. For polymer mixtures in the 
strong segregation limit, we have three length scales to consider: the length scale b of an 



effective bond, the radius of gyration R g = 5yiV/6, and the "intrinsic" interfacial width 
Wo, which is controlled by the Flory-Huggins parameter \- According to Helfand's [0, TA 
self-consistent field theory, this intrinsic width is 



(14) 



However, the finding of recent computer simulations that the \ parameter itself is 
not constant in the interfacial region in the strong segregation limit — due to subtle 
rearrangements in the chain configuration (enhancement of self- contacts, etc.) - has 
raised doubts on the accuracy of Eq.14, and in fact Eq.|l4| is not in agreement with the 

D . Thus, the proper estimation of w 



interfacial widths observed in the simulations 22 



and he 
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may be a problem. Semenov |2(| has suggested that wo should be used in q n 
concluded 

2 

Qmax • 

w 

Note the absence of a factor ir in Eq.|T3|: While in q min = 27r/L the constant 2ir simply is 
there due to the periodic boundary condition, there is no condition that actually would 
fix the constant of proportionality between q max and w ! 



3.3 Capillary wave broadening of the intrinsic profile: the con- 
volution approximation 

The treatment of the previous subsection has implied that the local interface is defined 
by the coordinate z int (x,y) = D/2 + h(x,y) only, i.e., the interface is considered as a 
sharp boundary without any intrinsic structure. While this "sharp kink" approximation 
is reasonable when long wavelength properties are concerned, it clearly does not make 
sense on small length scales of order u> : For length scales L of order 2ir/q max we would 
have q-max/flmin — 1; Eq.|13| predicts a vanishing mean square interfacial width, which 
is not a physically sensible result. The standard remedy || of this situation is to convolute 
the intrinsic profile p™^ (z) (which may be taken of the form of Eq|| in our case) with a 
Gaussian distribution of local interface heights, 
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such that the apparent profile becomes 
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For explicit calculations it is convenient to replace the tanh function by an error function 
with the same slope at the midpoint of the profile, 
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We now define the interfacial width w in terms of the maximum slope of the apparent 
profile, 



z=D/2 



Using Eq.18 in Eq.17 we find 
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(19) 
(20) 



and combining this result with Eq.[l3| yields 



w 




(21) 



Eq.|2l] is the well known result that the mean square broadening due to capillary waves 
simply has to be added to the square of the intrinsic width. 



3.4 Delocalized interfaces in confined geometry 



While Secs. p.2| , |Q| describe the situation of Fig. la in the limit D — > oo keeping L fixed, 
we now consider the inverse limit, keeping D fixed but letting L — > oo. Then the lower 
cut-off is no longer 2n/L, but rather a correlation length £\\ comes into play []3T], [S3). We 
briefly review the derivation of that length |33j. One again uses a Hamiltonian of the form 
of Eq.|iy, but now we need to explicitly consider potentials V(h) exerted by the walls on 
the fluctuating interface 
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or more explicitly 
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(23) 



Note that potentials of the form exp(— kD/2) cosh(Kh) simply arise from exponentially 
decaying forces due to the walls 

exp [— KZint] + exp [—k(D — Zi nt )} = 2 exp [— kD/2] cosh (ah) . 

Here the first potential changes its sign at a temperature T w to make allowance for the 
fact that in the semi- infinite system a wetting transition occurs at T w . The constant k in 
mean field theory |H| is identified with the inverse bulk correlation length while more 
refined treatments imply [[H], |60|, ^T], [62|, Q that k~ 1 gets renormalized by a factor (l+u/2) 
where us is the famous parameter entering the theory of critical wetting PBL |6l 



k- 1 =£(l + u /2), 



CO 



(24) 



Finally, a , bo are phenomeno logical constants that can be calculated from a microscopic 
theory of interfaces near walls for Ising models |]5S , |69| , [70| , but are not known for polymer 
mixtures, of course. 

Also the temperature dependence of the constant uo is known rather accurately for 
the Ising model ffl| , but not here. Fig. 2 shows a log- log plot of the expected behavior 



of uj vs N(x/Xcrit — !)• In the mean field critical regime of a polymer mixture (i.e. 



], |73| N\l — x/Xcrit\ 1) one predicts [§, |56| for a lattice 
model where every site is taken by either an A monomer or a B monomer (i.e, pb~ 3 = 1), 



I 1 - x/Xcrit\ < 1, but |56 
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and hence the mean field prediction for to is [f74|1 
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Very close to the critical point (i.e., for AT|1 — x/Xcrit\ <^ 1) one expects a crossover to 
the universal value of the Ising model fFI] ujj 



l stng 



0.86. However, more relevant for the 
present simulations is the behavior of uj in the strong segregation limit. While in the 
critical regime we have a simple relation between the (intrinsic) width of the interface wq 
and the correlation length [p6|l , wq = 2£, this is no longer true in the strong segregation 
limit (where £ ~ b\/~N /6 while w$ is independent of N, cf. Eq.[14]). Since in the critical 
region also could be expressed as uj = [iiWQcr} -1 , we suggest that this expression continues 
to be valid in the strong segregation limit rather than the expression given in Eq.|23], cf. 
Fig.2. Using pi, 111 a SSL 



b 2 Wx/6 and Eq.[14] one obtains 
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where in the last step the mean field result Xcrit = 2/AT was used. Consequently, we 
assume that also in the equation for k (Eq.|24D one should replace by 2/wo in the 
strong segregation limit. 

Expanding cosh(fi;/i) wl + (K,h) 2 /2 and omitting constant terms, Eq.p3| reduces to 



n, ff m = /^{|[(|) 2 +(|) 

exp (—kD/2) + 46 exp (-kD) 



T-T,„ 



n 2 h 2 



We see that the second term changes sign at [31], |33 | 



T C (D) = T u 



1 — 4 — exp 

a 



(2f 



(29) 



which is the critical temperature of the interface localization-delocalization transition [31 



33| . Of course, even on the mean field level the approximation Eq.^8] is only valid for 
T > T C (D). In this regime, Eq.^8] is formally identical to Ginzburg-Landau theories of 
second-order phase transitions 0, |3|, f|, and thus it is straight forward to obtain 
the correlation length £||, using the standard relation 



(30) 



to obtain 



.2 2a 2 / kD\ (T-T c (D) 



v=7^ -T hf" • (31) 



As expected, £|| shows a mean field type divergence at T = T C (D). Even for T T C (D), 
however, £|| is very large for high values of D , due to the factor exp(fc.D/4) in the relation 



£ll « k 1 yaT w /(2aoT)exp(K,D/4). If we now assume that it is the length £|| rather than 
L that cuts off the capillary wave spectrum, we conclude q min = 2tt/^\ and hence Eq.pl] 
becomes 

u> 2 = wl + — — h const, (32) 
16o~ 



the constant being const = (4a) -1 In n~ l q max ^ oT w j '(2a T) / '(2ir) . As discussed above, 
this constant is hard to estimate, and in the simplest approximation it has been ne- 
glected [0. This is justified in the mean field critical region, since ln[...] should be of 



order unity, and w 2 . 3> (4cr) _1 there, but it is not obvious that this approximation is 
accurate in the strong segregation limit. And again the question must be asked how reli- 
able it is to use ^\\ from Eqj3l| as a cut-off in a calculation with the "free" capillary wave 
Hamiltonian, Eq.[IT], rather than calculating s 2 directly from the full Hamiltonian, Eq.|2"5|, 
without further approximations. Also effects such as a possible position-dependence of 
the interfacial stiffness (the prefactor of the term (dh/dx) 2 + (dh/dy) 2 ) are neglected, 
although theories have suggested such effects for Ising-like systems ^8] |69], [7(J Hence a 
test of Eq.|32] by computer simulation is of significant interest. 



4 Interfacial profiles, interfacial widths, and the cor- 
relation length £|| 

Fig. 3 shows some raw data of order parameter profiles in a semi-grand-canonical simula- 
tion using L = 256 for films of various thicknesses D ranging from D = 16 to D = 64. 
The broadening of the profiles with increasing film thickness D is clearly recognized, and 



qualitatively these data are very similar to earlier results for the Ising model [j34| . From 
a fit of these data to tanh profiles, estimates for the apparent width w as function of L 
and D is obtained. 

In Fig. 4, this width w is plotted vs D for 32 < L < 512 in the semi-grand-canonical 
ensemble (sg) and the canonical ensemble with chain exchange (e) , respectively. We find 
three different regimes in our simulations. In the first one, D is small, and the data 
depend neither on L nor on the type of chosen ensemble. The increase of w 2 with D for 
D < 16 even is somewhat stronger than the expected (Eq.|32|) linear variation. This is 
due to the fact that for very small D (D < wq) not only capillary wave type fluctuations 
are suppressed, but even the intrinsic profile gets "squeezed" , and we expect w oc D (or 
w 2 oc D 2 ) as w — ► 0; there is just an essentially linear interfacial profile between such 
walls at extreme proximity of each other. In particular, for D — > the interfacial width 
w must vanish. Note that additional effects can also be expected due to the squeezing of 
the chains in quasi- two-dimensional configurations if D < R g . 



In the second regime for 24 < D < 56 and large enough L, data in the semi-grand- 
canonical ensemble do show the expected variation w 2 oc D + const. For small L, however, 
a crossover to the third regime with a steeper rise occurs. This crossover is clearly seen 
for L < 128. We expect that for L = 256 a similar upturn should occur for D = 64, but 
this is not observed due to the relative large statistical error of this last point (note that 
the exponential increase of the correlation length £|| with thickness D, Eq.|31~l, leads to 
a dramatic slowing down of all Monte-Carlo simulations, and this is responsible for the 
fact that a small increase of D causes a large increase of statistical errors). This upturn 
of the curves for large D and relatively small L is interpreted by noticing that we now 
enter a regime where L < £||, and then Eq.^2| is not valid: In a semi-grand-canonical 
simulation with periodic boundary conditions the dominating fluctuation is a mode with 
wavenumber q = 0, i.e., a uniform displacement of the interface as a whole ! Thus, there 
is a finite fraction of the thickness D where the interface can wander back and forth like 
a rigid straight interface, without feeling much effects of the walls. In the limiting case 
of D — > oo and L small, the interface is freely fluctuating, and thus Eq.|l9| gives a linear 
dependence w = D/2 pSfl . 

In the canonical ensemble, on the other hand, this third regime is very different. For 
L < £|| no significant further increase of w 2 with D occurs; rather w 2 then saturates at 
a plateau value because in the canonical ensemble a fluctuation of the mean interface 
position is suppressed, and the capillary wave spectrum now is again limited by L (q m i n = 



2tt/L) rather than Hence, instead of Eq.32 we now have from Eq.pT| 



w 2 = w 2 + _L i n L + _L i n ( SZOEL ) (33) 
4a 4a V 2vr J K J 

independent of D, which is the standard logarithmic variation of the interface squared 
width with the lateral linear dimension, as it has often been discussed in the literature |1|, 
0, |, |, H, §§. Considering the values of w 2 for D = 64 and L = 32, 64, and 128 
as independent from the finite thickness, we can extract from the slope of w 2 vs logL 
the interfacial tension a and find very good agreement with the simulations of Miiller et 



al ]2J, i.e., a = 0.015. 

In order to bring out more clearly these three regimes for the semi-grand-canonical 
ensemble, we compare in Fig. 5a the data for L = 256 (which cover the most extensive 
range of D and are least affected by the finite size effects associated with the smallness 



of L) with theoretical prediction, Eq.32. As pointed out in Sec. 3^, neither the choices 
of Wo nor the choice of k nor the choice of the additive constant is unique. But the 
choice made for the parameters in Fig. 5a — £ = 3.6 from Monte-Carlo simulations and 
Wo = Wq = 4.65 as determined in a self-consistent field calculation for our model [22 



- is at least plausible, even if the good agreement with the data for D > 40 is a mere 
lucky accident. In Fig.5b, a linear plot of w vs D checks for both linear regimes: In the 
first regime for small D, we find one straight line through the origin, including data where 
w < Wo and using w(D = 0) = 0. On the other hand, the crossover to the third regime 
for large D can be seen whenever £|| becomes comparable to L/2. Since ^\\ depends on D 
as £|| oc exp(nD/4) (Eq.[H]), it is clear that the location of the onset of the crossover shifts 
to larger D with increasing L while the first crossover (when w reaches its "intrinsic" 
value) is independent of L, of course. When L is very small (L = 32), the second regime 



(where Eq.|32] holds) is completely suppressed. For L > 64, this regime becomes more and 
more pronounced and an increasing number of data points fall on the expected L — > oo 
curve. Fig. 5 illustrates how difficult it is to give interfacial widths observed in simulations 
a precise meaning. 

We now discuss the correlation function g(r,z) (Eq.|j), see Figs. 6-8. In Fig. 6, an 
example for L = 256 and D = 32 in the semi-grand-canonical ensemble is shown. It is 
seen that near walls this quantity has a few oscillations for small r, reflecting the structure 
of the radial density distribution function, while in the center of the thin film (z = 16 in 
Fig. 6) a rather slow decay occurs. This slow decay is the hallmark of the "soft mode phase" 
described by the theory of Sec. |3.4| resulting from interfacial fluctuations. In Figs. 7 and 8, 
the correlation function g(r) = g(r, z = D/2) in the center of the film is studied. A 
reasonable fit to the expected form of the correlation function (Eq.[5[) is indeed possible. 
This is shown in Fig. 7 for different film thicknesses (16 < D < 48) using L = 512 and 
the semi-grand-canonical ensemble. The data points (symbols) are fitted to Eq.[5] (solid 
lines), and in this way the correlation length £|| is determined. However, we must add the 
caveat that correlation functions of the order of 10~ 3 and distances of the order of r > 40 
are clearly very hard to measure reliably in computer simulations. From Fig. 8 it is clear 
that simulations in the canonical ensemble underestimate g(r) at large r systematically 
if they are based on local moves and slithering snake moves only. Such observations as 
shown in Fig. 8 have led us to use the present canonical algorithm where also two identity 
exchanges of chains (A — > B and B — > A) are attempted as an additional move (data 
shown by squares in Fig.8). It is also disturbing, however, that there seem to be small 
systematic differences between the data for the canonical and the semi-grand-canonical 
ensemble (the estimates for £\\ from Fig.8 are £\\(sg) = 22.5 (circles) and £||(e) = 20.0 
(triangles), even though L = 256 exceeds C,\\ by about an order of magnitude. Clearly, 
for L —>■ oo both ensembles should yield identical results for £||, and in fact for L = 512 
and D < 40 this is found (cf. Fig.9). Presumably the observed discrepancy tells that the 
errors in £|| (due to statistical errors of g(r)) are of order of 10%, or there are unexpectedly 
large finite size corrections. 

Nevertheless it is gratifying that the data for £|| are compatible with the predicted expo- 
nential variation £|| oc exp(/tZ)/4), as shown in Fig.9. The same choice of parameters as as- 
sumed in Fig. 5a for the interfacial width (Eq.|32]) gives a good representation of £|| as func- 
tion of D. This, in turn, implies that a lot of the ambiguities about the choice of param- 
eters is avoided if we interpret Eqsj21| and |32| in the form w 2 = Wq + (4a) -1 \n(q max ^\\/27i) 
and simply use the measured £|| in this equation: a straight line with the same slope as 
included in Fig. 5a must result, irrespective of what one assumes for Wq, q ma x, and ui. 



5 A direct test of capillary wave concepts for inter- 
facial fluctuations 

The good agreement between the simulation results shown in Figs. 4, 5, 9 and theoretical 
concepts based on the capillary wave Hamiltonian suggests to study the interfacial fluc- 
tuations in more detail. In order to do this, it is convenient to split the system up into 
columns of block size B x B and length D and to determine in each column a local 



interface position h(x, y) using the concept of the Gibbs dividing surface (in a constant 
region in the center of the film in order not to be affected by bulk fluctuations). Fig. 10 
shows snapshot pictures of instantaneous configurations of the interface h(x,y) obtained 
in this way for a minimal coarse-graining B = 2 (columns have the width of the size of 
one monomer) and the physically more plausible choice B = 8 (comparable to the radius 
of gyration). While for B = 2 an extremely rugged "interfacial landscape" results, for 
B = 8 rather pronounced large-amplitude interfacial fluctuations are clearly visible. 

Having found h(x,y), we can proceed in the spirit of SecJOj and decompose h(x,y) 
into Fourier modes h(q), taking into account that possible q- vectors on the lattice are 
given as q = {q x ,q y ) = 2n/L with i,j integers. For simplicity, we consider only a 

wavevector in a lattice direction, and hence classify the h(q) as hi. Since Eq.[T2| implies 
(\hi\ 2 ) = (crq 2 )^ 1 oc i~ 2 , we present a log — log plot of (\hi\ 2 ) vs i~ 2 in Fig. 11. Pure 
capillary waves in this representation should show up as straight lines with a slope of 
unity. For L = 256 and D = 56 and 64, this is indeed found for 2 < i (< 20), but even 
for this largest studied thicknesses we still see deviations due to finite D for i = 1 ! For 
D < 48, this effect becomes more important and modes i — 2,3, ... are affected by the 
finite thickness D. One could suspect that there is a problem of equilibration of the largest 
length scales of interfacial fluctuations, but this is not a problem here as an analysis of 
the autocorrelation function Chhif) = [(hi(t)hi(0)) — (hi) 2 ]/[(h 2 ) — (hi) 2 } of the modes as 
a function of time measured in Monte-Carlo steps (MCS) shows (Fig. 12). The relaxation 
times of these fluctuations (r = O(10 A MCS)) are much smaller than the simulation runs 
(r = O(10 6 MCS)). 

An alternative analysis is to vary the block size B over a wide range and estimate the 
mean square width w 2 as a function of block size This is shown in Fig. 13 for L = 128. 
Here the interfacial position h is determined in all columns. The monomer profile is taken 
relative to this position, pa{z — h), and after averaging over all columns in the system, 
the interfacial width is determined by a fit to a tanh profile. For D = 56 and 64, one sees 
a constant slope for B > 8, as expected for capillary waves independent of film thickness 
D. Note that the largest block size B = 128 is different from the smaller ones. For B = L 
periodic boundary conditions apply. Thus, the capillary wave spectrum is cut off by 
Qmin = 2ir/L. In contrast, for B < L the boundary conditions are free (as in experimental 
situations) and a supplementary mode appears characterized by q = ir/B. For D = 16, 
there is an upper cut-off block size which gives a w independent of B, but typical for 
D — 16. For 24 < D < 48, we find a rather slow crossover from a D-dominated to a L- 
dominated interfacial width for the largest block size. On the other hand, the curves differ 
as well for all the small block sizes in function of the film thickness D, and this difference 
is almost constant for 4 < B < 32. This is a sign that there are two contributions to the 
variation of the interfacial width with B and D: The first arises from capillary waves; it 
depends on the block size B, and the finite thickness D acts as a new cut-off length scale 
£||. On the other hand, there is a second contribution that depends on D only, not on B. 
It can thus be attributed to a change in the "intrinsic width" of the interface due to the 
presence of the walls. 

For very small block sizes, w 2 becomes flat, and from this crossover one hence might 
attempt to estimate both w and q max directly. This would yield w ~ 3 and q max ~ 27r/5 
for D = 24, and wq ~ 4 and q m ax ~ 27r/6 for D = 64. But exact values are hard to 



give since there is a no sharp cut-off, but rather a smooth crossover. The three important 
length scales monomer size, correlation length £, and width w being close to each other, 
this analysis does not really clarify the questions concerning the cut-offs. 

Nevertheless the change of the intrinsic width with D extracted from Fig. 13 can 
be included quantitatively in our capillary wave analysis. For this purpose, we define 
w(D,B = 8) = wq(D) as the intrinsic width of our system. With this choice we are 
on the safe side, i.e. the further variation of w(D, B) in function of B can be described 
in terms of a capillary wave broadening. For B < 8, Fig. 13 indicates the onset of the 
crossover to a saturation value of wq. We take then the data from Fig.5b for w 2 (D) and 
in Fig. 14 w 2 (D) = w 2 (D) — Wq(D) is plotted. In that way, a constant slope of w 2 vs D 
is obtained over the whole range of simulations as required by Eq.S2[ The straight line 
shows the theoretical prediction with the same choice of parameter as before. The main 
result is that both slopes agree very well. Since in our case Wq CF ~ WQ m (D — > oo, B — 8), 
even the absolute values are the same (wq 1171 denotes the interfacial width extracted from 
simulations). Of course, from our data one could as well define w(D,B = 4) as intrinsic 
width. Then we would recover the theoretical result by choosing a constant term in Eq.[J^ 
larger than zero. 

Finally Fig. 15 considers the distribution of interface locations P(h). It is seen that 
Gaussian distributions give a very good fit to such data. Using w 2 (D,L = 128) = 
w 2 {D) + (n/2)s 2 (D) and w 2 {B = 8) from Fig.13, we find w 2 {D = 24, L = 128) = 
15.1 + (tt/2) 3.54 2 = 34.8 and w 2 (D = 64, L = 128) = 21.8 + (tt/2) 4.67 2 = 56.1, in good 
agreement with the direct results in Fig.13. This shows once more the self-consistency of 
our analysis. 



6 Conclusion 

In this paper we have studied the size-dependence of interfacial profiles between coexist- 
ing phases of mixtures confined between parallel walls, considering both the dependences 
on the linear dimension L parallel to the interface and the perpendicular linear dimen- 
sion D. The temperature of the system is above the wetting temperature T w , and the 
interface is delocalized and on average in the center of the film. We argue that in this 
regime the interfacial correlation length £|| oc exp(K.D/4), proposed first by Parry and 



Evans |!T] assuming that k is the correlation length £ of concentration fluctuations 
in the bulk, acts like a cut-off for the capillary wave spectrum, q min = 2vr /C| | - Us- 
ing the convolution of capillary wave broadening with the "intrinsic" profile of width 
w , we propose in the limit L ^> £|| the formula for the total interfacial squared width 
w 2 = Wq + (4a)" 1 ln(q ma x / <lmin) ~ Wq + (16cr) _1 kD , where a is the interfacial tension. 
This phenomenological approach is tested by a Monte-Carlo simulation of a model for a 
symmetrical polymer mixture in the strong segregation limit, and reasonable agreement 
is found, although the precise choice of parameters (cut-off q max , length capillary 
parameter ui in the relation k~ 1 = £(1 +u/2) proposed by recent theories |6l], [62|, |63| ) 



may need refinement. 

Of course, the above formula is only an asymptotic result for D ^> wq] for smaller 
values of D we find a regime where capillary wave type excitations are not yet important, 



and rather the intrinsic profile itself is "squeezed down" by the confinement. In this 
region w oc D (with a prefactor considerably smaller then unity) until w has reached 
roughly the value Wq of the intrinsic interface without confinement. Note that there D/2 
is comparable to the largest eigenvalue of the gyration tensor of the chains, and hence an 
effect on the chain configuration due to the confinement is already expected. This effect 
will be studied more closely in future work. 

The limit £|| > L deserves attention, too. We find that the above behavior w 2 ~ 
u>o + (16cr) _1 /t-D only holds until about £|| = L/2, and then a regime sets in which 
depends on the type of ensemble used: In the canonical ensemble, the average position of 
the interface is fixed at the middle of the system, and the capillary wave spectrum is cut 
off by q min = 2ix / L rather than q min = 27r/£||. As a result, w 2 = w 2 + (4a) -1 \n(q max L/27i) 
in this regime, independent of D. In the semi-grand-canonical ensemble, on the other 
hand, there is also a fluctuation described by a wavenumber q = 0, i.e., a uniform mode 
in which the interface fluctuates back and forth as a whole without bending. This leads 
to a behavior w cx D for D — > oo. 

As a result, the meaning of interfacial profiles and interfacial widths is rather delicate. 
This is not only a problem of computer simulations — it is encouraging to note that 
the analysis presented here is immediately useful for the interpretation of corresponding 
experiments |34|, |47|| . Of course, in experiments it may be sometimes necessary to work 
with long range surface forces (inverse power law decay) rather than the exponential decay 
assumed here. It is obvious that this will reduce the exponential growth of £|| with D into 
a power law growth as well, and since ln£|| = ln.D p = plnD for any power p, we obtain 
a behavior w 2 oc InD rather than w 2 oc D, as found here. A detailed analysis of this 
situation is still lacking. 

Many of our findings carry over to the other simulation geometries considered in 
Fig.l, such as antiperiodic boundary conditions or fully periodic boundary conditions. 
Note that in the latter case only the canonical ensemble is useful because otherwise the 
situation with two interfaces is unstable. In the antiperiodic case both ensembles may be 
used, but since one translates the interface to the center of the system by construction, 
the "uniform mode" leading to the above behavior w 2 oc D 2 is absent, and the finite 
size effects in the canonical and semi-grand-canonical ensembles are very similar then. 
We do hope that the present analysis will stimulate a careful assessment of these finite 
size effects in future simulations of interfacial properties since we feel that much previous 
work has been hampered by inadequate interpretation of fluctuation-broadened interfacial 
profiles as "intrinsic profiles" . The present work indicates that a unique characterization 
of "intrinsic profiles" encounters difficulties even if size effects are carefully taken care of. 
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Figure captions 



Fig.l: Typical geometries used in computer simulations (a,b) and in experiment (c) in 
order to study interfacial profiles. Only one lateral direction parallel to the average 
interface plane is shown. The schematic plots indicate "coarse-grained snapshots" 
of the system, i.e., bulk fluctuations in both the A-rich phase and the B-rich phase 
are averaged out, and the resolution perpendicular to the interface does not resolve 
the "intrinsic" interfacial profile in these pictures. Case (a) shows geometries where 



a single interface occurs |23], [3^, |33|, |34], [39|], while case (b) refers to a situation with 



two (on average parallel) interfaces confining a slab like region of B-rich phase in 



between [ |0| . While the left case of (a) has the disadvantage that perturbations due 
to the two hard walls must affect the interfacial profile if D is not large enough, 
the right case of (a) is only possible for strictly symmetric mixtures ( "antiperiodic 
boundary conditions" means that an A-particle that leaves the simulation box on 
the right side reenters as a B-particle on the left side, while no identity switches 
of particles occur across periodic boundaries). In case (c) it is assumed that the 
direction of gravity is from left to right and that the B-rich phase is favored by the 
air while the A-rich phase is favored by the right wall (walls of the container at the 
bottom and at the top are assumed to be neutral). Note that here linear dimensions 
are not to scale (e.g. the sample lateral dimension L samp i e may be macroscopically 
large and exceed the lateral resolution length scale L by several orders of magni- 
tude), and also the fluctuations of the coarse-grained interface are simplified (short 
wavelengths being omitted, amplitudes of long wavelengths exaggerated for clarity). 

Fig. 2: Log-log plot of u vs x/Xcrit — 1, for chain lengths ranging from N = 128 to 
iV = 1024, showing the predictions for ui S i ng , ujmf (Eqj26|), and ujssl (Eq.^7]) us i n g 
pb 3 = 1. 

Fig. 3: Order parameter profile m(z) vs z — D/2 for films of thicknesses D = 16,32,48, 
and 64, as indicated by different symbols. The interface is centered at the origin 
z = D/2, walls being situated at z = and z = D. All data refer to a lateral system 
size of L = 256 and the semi-grand-canonical ensemble. The lines through the data 
points are fits to m(z) = mfetanh \{z — D/2)/w]. Note that statistical errors are 
smaller than the size of the symbols. 

Fig. 4: Plot of the square of the apparent interfacial width w 2 vs film thickness D for 
several choices of L (L = 32, 64, 128, 256, and 512) and two choices of statistical 
ensembles: the canonical ensemble including identity exchanges (opaque symbols 
marked by (e)) and the semi-grand-canonical ensemble (filled symbols marked by 
(sg)). Note that the accuracy of these data is very good for small D, while for the 
two largest values of D (D = 56, 64) statistical errors are estimated to exceed the 
size of the symbols several times. For L = 512, values of w 2 are determined up to 
D = 48 and, within statistical error, no difference between the ensembles is found 
and, thus, points in Fig. 4 lie on top of each other. 

Fig. 5: a) Squared interfacial width w 2 vs film thickness D for L = 256 in the semi-grand- 
canonical ensemble. The straight line shows Eq.p2l neglecting the last term on the 



right hand size and using k 1 = £(1 + ui/2), £ = 3.6, uj 1 = 47r£ 2 a, a = 0.015, and 
Wo = w scf = 4 65 n _ 

b) Interfacial width w plotted vs D for L = 32, 64, 128, and L = 256 using the data 
from the semi-grand-canonical ensemble simulations. For small D, one has a linear 
variation w oc D up to w ~ u7o (squeezed intrinsic interface). Arrows show that 
the regime of free interface fluctuation sets in for £|| = L/2. The dashed line is the 
limiting case of completely free fluctuations of the interface, i.e., w = D/2. Also 
the self-consistent field estimate for the intrinsic interfacial width (wq CF ) is shown. 
Data for £n have been taken from Fig. 9 below. 

Fig. 6: Correlation function g(r,z) (Eq.|j) plotted as function of both variables r = 
\/x 2 + y 2 and z for the case L = 256, D = 32. 

Fig. 7: Semilog plot of correlation functions g(r) = g(r,z = D/2) in the center of the 
thin film vs r for L = 512 and general choices of D as indicated. All data refer 
to the semi-grand-canonical ensemble. Curves are fits to g(r) to the form g(r) oc 
r -1 / 2 exp(— r/£||), cf. Eq|| 

Fig. 8: Semilog plot of correlation functions g(r) vs r for L = 256 and D = 32 for different 
ensembles: semi-grand-canonical ensemble (circles), canonical ensemble (triangles) 
including chain identity exchanges, while squares indicate results with local moves 
and slithering snake moves only. Curves are fits to Eq.|5|. 

Fig. 9: Correlation length £n plotted vs D obtained by a fit of g(r) to Eq.[| for L = 512 
and film thicknesses D < 48 using both data for the canonical ensemble (triangles) 
and the semi-grand-canonical ensemble (circles). Within statistical errors, both 
ensembles yield identical results for D < 40, but differ for D = 48. Full curve is the 
formula £|| = £ exp(/?.D/4), using k" 1 = £(1 + uj/2) and £ = 3.6 as obtained from an 
independent simulation. 

Fig. 10: Typical plot of a snapshot picture of the local interface positions h(x,y) for 
D = 64, L = 64 using the minimal coarse-graining (i.e., block size B — 2) in case 
(a) and modest coarse-graining (B = 8 ~ R g ) in case (b). 

Fig.ll: Log-log plot of the mean square Fourier component (\hi\ 2 ) of the interfacial 
fluctuations vs i~ 2 for L = 256 and the canonical ensemble. Several choices of D are 
included as indicated in the figure. The straight line shows the theoretical prediction 
from the capillary wave Hamiltonian, i.e., (\hi\ 2 ) = (aq 2 ) 1 = [L / (2%i)] 2 a^ 1 . 

Fig. 12: Semilog plot of the autocorrelation function Chhif) = [{hi{t)hi{Q)) — {hi) 2 ]/[{h 2 ) — 
(hi) 2 ] vs time t (in Monte-Carlo steps MCS). Data refer to L = 256, D = 48 
using the canonical ensemble with chain identity exchanges. The two slowest modes 
(i = 1,2) are shown. Lines are fits of — ln[C^(t)] to an exponential decay. Even for 
the slowest mode (i = 1) the autocorrelation time r « 1.5 10 4 MCS is much smaller 
than the time of typical simulation runs. 



Fig. 13: Plot of the squared interfacial width w 2 vs the block size B (note the logarith- 
mic scale) for a system of lateral linear dimension L = 128 and thicknesses D as 
indicated. 

Fig. 14: Same date as presented in Fig.5a (L = 256, semi-grand-canonical ensemble), but 
only the contribution w 2 (D) = w 2 (D) — Wq(D) due to capillary waves is shown. 
The intrinsic width Wq(D) is taken from Fig. 13 (L = 128, canonical ensemble) as 
the width determined in blocks of B = 8 ~ R g . The straight line is the broadening 
according to theory (cf. Eq.p2|), i.e., w 2 = kD/(16ct). 

Fig. 15: Distribution of interface locations P(h) for L = 128, 24 < D < 64 using a coarse- 
graining size B = 8. Lines are fits to the theoretical assumption of a Gaussian 
distribution P(h) = (2txs 2 )~ 1 I 2 exp[— ti 2 /(2s 2 )]. The values of s obtained by this fit 
are indicated. 
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